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Abstract. We report here the discovery of a 30"-chain of embedded Herbig-Haro (HH) objects in the p Ophiuchi 
dark cloud. These HH objects were first detected during a deep Ks-hand observation (completeness magnitude 
for point source~19) made with NTT/SOFI. We confirm their nature with follow-up observations made with 
H2 11=1-0 S(l) narrow-band filter. We argue that they belong to two different jets emanating from two Class I 
protostars: the main component of the recently resolved subarcsecond radio binary YLW15 (also called IRS43), 
and IRS54. We propose also to identify the [S H] knot HH224NW1 (Gomez et al. 1998) as emanating from a 
counterjet of YLW15. The alignment between these HH objects and the thermal jet candidate found in YLW15 
by Girart et al. (2000) implies that this jet is not precessing at least on timescale ~(2-4) x 10'* yr. 
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1. Introduction 

In the 1950's, objective prism surveys of dark clouds re- 
vealed optical small nebulae, associated with young stars, 
showing emission line spectra with very weak continua. 
Today, these Herbig-Haro (HH) objects are interpreted as 
shocks produced by the interaction of outflows from Young 
Stellar Objects (YSOs) and the interstellar medium (see 
review by [Reipurth fc Raga 1999| ). 

The p Ophiuchi dark cloud is one of the nearest 
(d~140pc) active site of low-mass star formation, dis- 
playing a rich embedded cluster of ~200 YSOs (see the 



ulae a relatively high [Sll]/Ha ratio, which is a tracer of 
low excitation conditions, and thus indicates that H2 is 
not destroyed after the bow shock; they suggested to use 
the shocked-H2 transition at A=2.121/xm to probe deeply 
embedded HH objects. Up to now in the p Ophiuchi dark 
cloud, shocked-H2 imaging was only performed to study 
molecular shocks in the CO out flow of the Class pro - 
tostar prototype VLA1623 (e.g., [Davis fc Eisloeffel 1995|), 



upda ted census by the ISOCAM survey; Bontemps et al 
2001 ) . It is thus one of the most suitable locations for the 
search of HH objects. Optical emission-line surveys (Ha, 
[SII]) have led to the detection of only 10 bona fide HH 
objects dReipurth fc Graham 1988|; [Wilking et al. 1997|; 



Gomez et al. 1998; Reipurth 1999). As outlined by Wilking 



et al. and Gomez et al., most of these HH objects are lo- 
cated at the periphery of the cloud, in the lowest extinc- 
tion area: optical surveys do not have access to embedded 



HH objects. Wilking et al. (1997) noted for all these neb- 
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which led to the detection of several H2-knots; one of th em 
was later detected in [SII] (HH313; Gomez et al. 1998 ). 

We report here deep near-IR observations of the p 
Ophiuchi dark cloud unveiling bow-shape structures and 
knots, and complementary observation with narrow-band 
filter confirming their nature as H2 shocks. As these ob- 
jects are not optically visible, we called them embedded 
HH objects, and discuss the possible exciting sources. 



2. Observations and results 

During the period April 4-7 2001, deep J, H, ifg-band ob- 
servations of the p Ophiuchi dark cloud were made by N.G. 
with NTT/SOFI, as follow-up of Chandra and XMM- 
Newton X-ray observations. The complete report of this 
follow-up, consisting of 5 pointings, will be published in an 
upcoming paper (Grosso et al., in preparation). We will 
focus here only on the pointing related to the detection of 
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Fig. 1. i^s-band finding chart based on 2MASS survey im- 
age, showing the NTT/SOFI pointing field of view. The YSOs 
of this area are marked according to their IR classifications 
(diamond/square/cross for resp. Class I/II/III IR source; see 
Lada 1991) from the ISOCAM mid-IR survey (Bontemps et al. 
2001). The box indicates the area enlarged in Fig. 

new embedded HH objects , and located SE of the milli- 
metric dense core Oph-B2 ( Motte et al. 1998 ; see Fig. |l|). 
This area was already surveyed by optical emission-line 
surveys, but no HH objects were found. 

We used the auto jitter imaging mode to obtain a to- 



of June 8 2001 (Fig. ^.b). The flat-field was computed from 
the median of the dark substracted frames. The good see- 
ing conditions (FWHM~0.7") of this shorter observation 
gives much more details than the ifg-band image. 

To differentiate between pure line emission knots and 
continuu m emission from stars (inclu ding scattered light) 
we follow Davis fc Eislocffel (1995) who noted that the 
narrow-band filter reduces the intensity of the field stars 
by the ratio of the filter bandpasses (~10), whereas the 
intensity of the H2 knots is only reduced by ~2 (indeed 
spectra of H2 molec ular shocks have other emission lines 
in the iCs-band; see Smith 1995 ). We scale down the Ks- 
band image by the ratio of the filter bandpasses and the 
ratio of the seeing to have continuum features appearing 
with the same brightness in Fig. ^ and ||b. On the other 
hand, H2-shocks appear ~5 times brighter in Fig. ^3 than 
in Fig. ||a and can thus be easily identified. The H2-knots 
detected are labelled on the H2 contour map (Fig. ||c). 
Table |l| gives the position of these knots with photometry 
for a l"-aperture from the 5CT-filtered images. Before ap- 
plying wavelet filtering the H2 image was convolved with a 
Gaussian filter to have the same seeing as in the i^^s-band 
image. To obtain the H2 magnitudes, we scaled up the H2 
intensities by the filter-band width ratio and applied the 
magnitude zero-point of the Ks image. 

The knots A display a bow-shape structure, character- 
istic of low-speed bow-shocks which produce such arcs and 



limb- brightened structures of conical appearance (Smith 



tal integration time of 40 mn, with 6 s detector integration 1991). We also note in the H2 image some diffuse emission 



time (DIT) and 1 mn-exposure frames. The near-IR stan 
dard AS29-1 ( Himt et al. 1998| ) was observed several times 
during the night. We took five 6s-darks at the end of the 
night. The different stages of the usual IR data reduction, 
namely sky estimation and subtraction, frame recentering 
and stacking, including also dark substraction and flat- 
field division (we used the flat-field provided by the NTT 
team for our observational period), were performed using 
the j itter routine of the ESO's eclipse package (version 
3.8-1). Zero order astromc tric correction wa s applied us- 
ing the 2MASS field stars ( |Cutri et al. 2000|) as reference 
frame, reducing the offset residual to 0.7". The magnitude 
zero-point was obtained from our calibration star. 

Fig. Ija is an enlargement of the observation ob- 
tained with the Ks broad-band filter (A=2.162/zm, 
AA=0.275/im), unveiling a complex object showing bow- 
shape structures and possible emission knots. To detect 
weak structures in the J and H-hand images, we filtered 
it with the wavelet transform (MR/1 package; 3tarck el 
al. 199^ to a 5a significance threshold. In this manner, a 
few of knot candidates are also detected in the 7?-band, 
and only one is detected in the J-band. The shape of this 
complex object is reminiscent of the ones found in HH 
objects. If this HH object interpretation is correct, this 
iiTs-band detection must be mainly due to the 2.121 /xm- 
line of the H2 shocked emission. To confirm it, J. A. made 
a 20 mn-exposure observation (DIT=30s, 1 mn-exposure 
frames) with NTT/SOFI using the H2 v=l-0 S(l) narrow- 
band filter (A=2.124/im, AA=0.028^m) during the night 



downstream of the bow-shock which may be related to 
the excitation of H2 in the pre-shock region of a J-shock. 
We thus propose the knots A as the leading part of a jet 
coming from the SW direction, and outlined by the up- 
stream knots B. By contrast the knots C are displaced 
from the direction of this jet, and moreover the knots Ci- 
C4 are clearly elongated, roughly along a NS direction. We 
propose to explain these features by shocks coming from 
the East and producing the arc C2-C1-C5, with the down- 
stream knots C3 and C4. We will discuss in the following 
section the possible exciting sources of these two jets. 

To obtain an estimate of the H2-shock velocity we com- 
pare the observed near-IR photometry with the predic- 
tions for planar molecular shocks ( ^mith 1995| ). We up- 
date first these results for the NTT set of filters. We ex- 
tract the filter transmission profiles from the plots pro- 
vided on the SOFI web page, and convolve them to the 
unfiltered molecular shock line fluxes (Michael D. Smith, 
private communication). Zero points ior J — H and H — Ks 
colors are derived by requiring the colors of AO spec- 



trum template (Pickles 1998) to be zero. We observe 
f(irs)/f(H2)~1.4 and -2.0 for the leading part of the 
two jets, which must be compared to f(ii'g)/f(H2)=2.0- 
3.3 {resp. 1.9-3.3) for J-shocks {resp. C-shocks) with ve- 
locity range 8-22kms^^ {resp. 20-45kms^^). This im- 
plies J-shock {resp. C-shock) velocity <10kms~^ {resp. 
20kms~^), well below the H2 dissociation velocity in 
molecular cloud (^^22, — 47kms^^ for resp. J, C-shock; see 
[Smith 1995| ). Hence, the H2-hne emission is really map- 
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Fig. 2. Deep near-IR images of the embedded HH objects discovered with NTT/SOFI. (a) ATs-band image scaled so that 
continuum features from stars appear as bright as in the H2 !;=l-0 S(l) narrow-band filter image. The pixel size is 0'.'29, and the 
grey color scale is linear. J-band {resp. iif-band) detections above 5a levels are shown with dashed {resp. continuous) contour 
levels (24.5, 23.5, 23.0 and 21.8, 21.2, 20.9 mag arcsec"^). (b) H2 t)=l-0 S(l) narrow-band filter image (including the continuum 
emission). H2 line emission features appear brighter in this narrow-band image than in the previous broad-band one. (c) H2 
t;=l-0 S(l) contour map (18.3, 17.6, 17.2, 16.9, 16.6, 16.5, 16.3, 16.1, 16.0, 15.9 magarcsec"^). The H2-knots are labeled by 
decreasing H2 magnitudes (see Table |l|). The first {resp. second) arrow joins the barycenter of the knots A-B {resp. C) and the 
Class I protostar YLW15 {resp. IRS54). The first arrow corresponds also to the symmetry axis of the knots A-B. 



ping the bow-shock; we can exclude the scenario where 
this emission would come only from the low-velocity wing 
of a fast dissociative bow-shock located downstream, with 
strong iron Hues in the J and H bands tracing the apex. 
These low-velocity embedded H2-shocks are reminiscent 



Table 1. H2-knot positions and magnitudes. 



of the low excitation HH objects already observed in the 
optical in this dark cloud. From the intrinsic color of low- 
velocity H2-shocks, ( J — iJ)n^0.8, and as suming the red- 
dening law quoted by Cohen et al. (1981) , we estimate the 
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extinction for knot Ci: Av=9.09x[(J-iJ)-(J-if)o]~23. 
This large visual extinction explains the non-detection by 
previous optical surveys. It will be possible to measure di- 
rectly Ay by J, i7-band spectroscopy using [Fell] lines at 
1.25 iiui, 1.64 /xm, which arise from the same atomic upper 
level. 



3. Discussion on the exciting sources 

The optical survey of HH objects in Taurus has shown 
that their frequency decreases rapidly with the age of the 



YSO, from Class I sources to Class H sources (Gomez et 
|al. 1997 ), i.e. from evolved protostars to classical T Tauri 
stars. This result is consistent with the decrease of the CO 
outflows observed from the C lass sources to the Class I 
sources ( Bontemps ct al. 199(i| ). To find the exciting source 
candidates of the HH objects reported here, it is thus rea- 
sonable to look for YSOs with IR excesses. On the basis 
of the HH feature morphology, we proposed in the previ- 
ous section shocks coming from the East {resp. SW) to 
explain the shape of knots C {resp. A-B). We checked by 
constructing a color-color diagram of the sources detected 
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in our deep observation, that there is no new embebbed 
YSO with IR excess in these directions. Two known YSOs 
are Eastward (see Fig. the Class II source GY350, and 
the Class I protostar IRS54. We propose to associate knots 
C with IRS54, the agreement with the knots C shape look- 
ing better (see Fig. Qc). Strong II2 v=l- S(l) line cmis- 
sion was detected in the IRS54 spectrum ( Greene fc Lada 
is unresolved in our H2 image, which shows only 
scattered light Eastward. To our knowledge this source 
has never been included in a CO outflow survey. 

Fig. H.c shows that the Class I protostar YLW15, also 
called IRS43, is on-axis of the knots A-B, ~10' away 
(0.4 pc for d=:140pc). This evolved protostar is a radio 
source, and was recently announced to be a subarcsecond 
radio binary ( Girart et al. 2000 ). The main radio compo- 
nent, YLW15-VLA1, is spatially extended with a position 
angle of 25°±2°(one sigma error), and Girart et al. pro- 
posed it to be a thermal radio jet candidate. The position 
angle (PA) of the knot chain, 22.3°, is compatible with 
the PA of this jet candidate within errors (Fig. |). As the 
probability for a coincidence by chance between these PAs 
is very low (2xl.645x2°/180°«0.037), the identification 
of YLW15 as the exciting source of knots A-B is likely. 
Fig. H displays the position of the HH objects of this dark 
cloud. The [SII] knot HH224NW1 ( pomez et al. 1998D 
was associated with the complex HH object HH224. We 
note however that HH224NW1 is displaced from the axis 
of HH224, and that the PA of YLW15 is 23.4°, thus also 
compatible with the orientation of the thermal jet candi- 



mal jet in this object. Moreover the alignment with the HH 
objects moving at V <10-20kms~^ implies that this jet 
is not precessing at least on timescale ~(2-4)xl0^yr. 

These observations confirm that HH objects excited by 
Class I protostars may be hidden by large extinctions in 
the optical, but can be easily unveiled and studied by deep 
near-IR observations. Embedded HH objects are probably 
very common, and it is essential now to obtain a reliable 
census, for instance to quantify the role of outflows in 
maintaining a high level of tur bulence in molecular cloud s 
and regulating star formation ( Matzner fc McKee 2000| ). 
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